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ABSTRACT 

Aims. We have analyzed high-resolution and high signal-to-noise ratio optical spectra of nearby FGK stars with and without detected 
giant planets in order to homogeneously measure their photospheric parameters, mass, age, and the abundances of volatile (C, N, and 
O) and refractory (Na, Mg, Si, Ca, Ti, V, Mn, Fe, Ni, Cu, and Ba) elements. Our sample contains 309 stars from the solar neighborhood 
(up to the distance of 100 pc), out of which 140 are dwarfs, 29 are subgiants, and 140 are giants. 

Methods. The photospheric parameters are derived from the equivalent widths of Fei and Feii lines. Masses and ages come from 
the interpolation in evolutionary tracks and isochrones on the HR diagram. The abundance determination is based on the equivalent 
widths of selected atomic lines of the refractory elements and on the spectral synthesis of C 2 , CN, Ci, Oi, and Nai features. We apply 
a set of statistical methods to analyze the abundances derived for the three subsamples. 

Results. Our results show that: i) giant stars systematically exhibit underabundance in [C/Fe] and overabundance in [N/Fe] and 
[Na/Fe] in comparison with dwarfs, a result that is normally attributed to evolution-induced mixing processes in the envelope of 
evolved stars; ii) for solar analogs only, the abundance trends with the condensation temperature of the elements are correlated with 
age and anticorrelated with the surface gravity, which is in agreement with recent studies; Hi) as in the case of [Fe/H], dwarf stars with 
giant planets are systematically enriched in [X/H] for all the analyzed elements, except for O and Ba (the former due to limitations 
of statistics), confirming previous findings in the literature that not only iron has an important relation with the planetary formation; 
and iv) giant planet hosts are also significantly overabundant for the same metallicity when the elements from Mg to Cu are combined 
together. 
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1. Introduction 

The fact that stars hosting a giant planet are, on average, 
more abundant in iron than stars in the solar neighborhood 
for which no planet has been detected (ISantos et ^ l200ll 
I 2 OO 4 I: iFischer & Valentill2005l: lGonzalezll2006l) is well accepted. 
Actually, it has been shown that this behavior is not exclusive to 
iron but it is sha r ed by several other metals (iBond et ai]l2006t 
iGilli et all I2006I: iNeves et aH I2009I; lAdibekvan et al.l l2012bl) . 
Other studies also suggested that this kind of anomaly may 
not only involve the metal content of heavy elements but also 
some light elements such as carbon, nitrogen, an d oxygen 
(lEcuvillon et al .1120041 l2006t iPetigura & Marc^Ol ill . 

Some authors have also been searching for differences in the 
abundance trends with metallicity for several elements in stars 
with and without planets. In other words, these works investi¬ 
gated whether the planetary formation in some way affects the 
observed abundances by compar ing stars in the same metallic¬ 
ity range. iRobinson et al.l (l2006h reported that the [Si/Ee] and 
[Ni/Ee] abundance ratios in planet-host stars are systematically 


enhanced over their comp arison sample of stars witho ut plan¬ 
ets of the same metallicity. iDelgado Mena et ^ (120101) derived 
the opposite behavior for Mg, and they found no differences re¬ 
garding [C/Ee], [O/Ee] and [Si/Ee]. In any case, they stated that 
their res ult for Mg d i sappea rs when only solar analogs are con¬ 
sidered. iNeves et^ (l2009l) also reported no significant differ- 
enc es in [X/Ee] for 12 eleme nts in stars with and without plan¬ 
ets. iGonzalez & LawsI (l2007l) derived lower [Al/Ee] and [Si/Ee] 
but highe r [Ti/Ee] in stars with p lanets in the range of metal- 
rich stars. iBrugamver et alJ (1201 lb . on the other hand, confirmed 
the [Si/Ee] enhancements in planet hosts, in spite of no [O/Ee] 
enhancemen ts. Overabu n dance s of s everal elements were als o 
reported by iKang et ^ (1201 ll) and lAdibekvan et alJ (1201 2ah . 
but only in t he range of metal- poor stars, except for Mn. Eor 
this element, iKang et ^ (1201 ll) found an enhancement in the 
whole range of metallic ity, a result that is not supported by 
lAdibekvan et alJ (l2012al) . The situation is, therefore, vast and 
varied, and also controversial for some elements. 

In a previous paper (Ida Silva et al.ll201ll from now on re¬ 
ferred to as Paper 1) we determined photospheric parameters and 
carbon abundances for a sample of 172 dwarf, subgiant, and gi¬ 
ant stars (out of which 18 with planets) using sp ectra available 
in the ELODIE database (iMoultaka et al.lU004t) . Our previous 
results did not point out, for instance, any significant difference 
between the carbon abundances of stars with and without plan¬ 
ets, for the three substellar groups. However, the analyzed sam- 


Send offprint requests to: R. da Silva, 
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* Based on public data from the ELODIE archive (Moultaka et al. 
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** Tables 4-8 are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 
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Fig. 1. HR (left panel) and color-distance (middle panel) diagrams for our sample of 140 dwarfs (□), 29 subgiants (A), and 140 
giants (O). The derived stellar parameters are also shown in a log,g vs. Tefr diagram (right panel) for comparison. Filled symbols 
represent stars with detected planets. The distance limit of 100 pc is represented by the dotted line. 


pie was quite small to provide conclusive results, and stimulated 
the inclusion of additional stars and other chemical elements. 

In the current work we have 309 FGK stars (out of which 
47 with detected planets) homogeneously analyzed in order to 
measure their photospheric parameters, masses, ages, and the 
abundances of volatile (C, N, and O) and refractory (Na, Mg, 
Si, Ca, Ti, V, Mn, Fe, Ni, Cu, and Ba) elements. The sample in¬ 
cludes dwarf, subgiant, and giant stars, mostly belonging to the 
Galactic thin disc. We have applied a set of statistical methods 
to the three subsamples, exploring their multivariate properties 
and, in particular, uncovering significant differences between 
dwarf stars with and without a known giant planet. In addition, 
we have searched for relations involving the stellar parameters, 
the abundances, and the condensation temperature of the ele¬ 
ments. This approach has been very often used in order to check 
whether any difference in the amount of refractory elements with 
respect to volatile ones has an influence on the formation of plan- 


etarv systems (see. e.ff.. Smith et al.l200lHMele 

ndez et al.l2009l 

iRamirez et al.ll2009HGonzalez Hernandez et al.l 

201 011201 3aybh. 


In Sect. |2] we describe the observations and the procedure of 
data reduction. In Sect.|3]we outline the methods used to derive 
the photospheric parameters and the chemical abundances. The 
determination of the evolutionary and the kinematic parameters 
is presented in Sect. |4] The results are presented and discussed 
in Sect. |5] and in Sect. in which we describe the statistical 
analysis. Finally, we present our final remarks and conclusions 
in Sect.|2] 


2. Observation data and reduction 

In this work, as a complement to our previous analysis of 172 
stars published in Paper 1, we have enlarged our stellar sample. 
Now we have 309 FGK stars from the solar neighborhood that 
include 140 dwarfs (out of which 31 with detected planets), 29 
subgiants (4 with planets), and 140 giants (12 with planets). 

The stars were observed with th e ELODIE high-resolution 
spectrograph (iBaranne et al.l 1 19961) of the Haute Provence 
Observatory (Prance), and the data are publicly available in the 
online ELODIE database. The ELODIE spectrograph provides 
spectra with resolution R ~ 42 000 in the wavelength range 
3895-6815 A. The criteria of selection are similar, but not equal, 
to those applied before since we have now considered a smaller 


inferior limit for the signal-to-noise ratio (S/N). Pollowing this 
new criteria we selected; 

i) stars for which the averaged spectra have S/N > 150; only 
individual spectra with a S/N > 20 and with an image type 
classified as object fibre only (OBJO) were used; 

ii) stars within a distance < 100 pc (parallax tt > 10 mas) and 
with spectral type between P8 and Ml; 

Hi) stars having no close-in binary companion that could con¬ 
taminate the observed spectra; we used the information of 
the angular separation betw een compon ents (rho) available 
in the Hipparcos catalogue (IESAII19^ and we chose only 
the systems with rho >10 arcsec; additionally, we searched 
in literature for any further information that could unveil a 
contaminant star in the observed field; 

iv) stars with {B - V) measured by Hipparcos and with cross¬ 
correlation parameters available in the ELODIE database; 
both (B - V) and the width of the spectral cross-correlation 
function are required to carry out a first estimate of the pro¬ 
jected rotation velocity vsin i of each star; and 

v) stars for which the determination of the photospheric param¬ 
eters is reliable (see Sect. 13.11) and that passed the quality 
control of the spectral synthesis (see Sect. 13.2b . 

A total of 1485 spectra of 309 stars has been analyzed. The 
subsamples of dwarfs, subgiants, and giants are plotted in the 
HR and color-distance diagrams of Eig. [T] As in Paper 1, we 
chose to classify as subgiants those stars situated 1.5 mag above 
the lower limit of the main-sequence and having > 2.0 mag. 
Note that the distance of dwarfs and subgiants is not limited to 
100 pc, but to about 80 pc, which represents a selection effect of 
the ELODIE observation surveys. Eor comparison, in the right 
panel of this figure we plot the logg vs. Tgs diagram with the 
stellar parameters derived in the current work. 

The spectra were reduced using IRAEQ routines for order 
identification and extraction, background subtraction, flat-field 
correction, wavelength calibration, radial-velocity shift correc¬ 
tion, and flux normalization (a global normalization first, and 
then a more careful one performed around regions containing 
molecular bands and atomic lines used in the spectral synthesis). 


* Image Reduction and Analysis Facility, distributed by the National 
Optical Astronomy Observatories (NOAO), USA. 
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Fig. 2. Synthesis of three spectral regions: C 2 (/15165, left panels) and CN (/i4215, middle panels) molecular band heads, and the 
O atomic line at 46300.30 (right panels). Two dwarf stars with different values of spectral S/N are shown: HD 168603 (S/N = 156, 
Teff = 5399 K, [Fe/H] = 0.11 dex, top panels) and HD 190771 (S/N = 299, T^ff = 5817 K, [Fe/H] = 0.14 dex, bottom panels). These 
are the S/N per pixel computed at 45550. Six spectra computed for different abundances (the N abundance is changed in the case of 
the CN band) and separated by 0.1 dex are plotted. The differences between observed and computed spectra (O-C) are also plotted 
in the bottom of each panel. 


3. Spectroscopic analysis 

3.1. Photospheric parameters and equivalent widths 


In this work, as in Paper 1, we performed a differential analysis 
relative to the Sun. The method requires a model atmosphere 
for each star, which in turn is computed using the so-called 
photospheric parameters: effective temperature Tsg, metallicity 
[Fe/H], surface gravity logg, and microturbulent velocity f. We 
computed these parameters through the excitation equilibrium 
of neutral iron and the ionization equilibrium between Fe i and 
Fe II lines, using model at mospheres derived by interpolation in 
the grid of lKuruc3 ( 1993^^s before, we used the MOOG LTE 
radiative code ( Snedenl 2002 h and the same routine developed by 
us to perform an automated calculation of the parameters and 
their uncertainties, allowing a fast and homogeneous approach. 


To compute the solar model atmosphere we adopted Tetr = 
5777 K, logg - 4.44, 4^ - 1.0 km s~\ and the standard solar 
abundances of ! Anders & Grevessel(ll989h . but with log e 0 (Fe) = 
7.47 to be consistent with the value adopted in Paper 1. We also 
assumed the Unsold approximation multiplied by 6.3 to account 
for the van der Waals line damping. This is useful for a better 
treatment of strong lines close to the spectral windows used in 
our spectral synthesis analysis. 


The atomic line parameters (wavelength, statistical weight 
of the lower level multiplied by the oscillator strength gf, and 
lower-level excitation potential ;t') of 71 Fe i and 12 Fe ii lines are 
listed in Table |2l The table also lists the parameters of the other 
elements analyzed by means of their equivalent width (EW). 
Table [3] contains the gf values of the elements with important 
hyperfine structure (HFS). These parameters were i nitially taken 
from the Vi e nna Atomic Line Database - VALD ( Ku pku_et^ 
I 2 OOOI 1 19991: iPiskunov et al.lll99^ iRvabchikova et alJll997h in 
the case of the elements without HFS, and from SteffenI (Il985h 
for those with important HFS. The gf values were then revised 
in order to fit the equival ent widths measure d in the spectrum 
of the Solar Flux Atlas (iKurucz et al.l Il984l) degraded to the 
ELODIE resolution. Eor the Sun and for all the other stars of our 
sample, the EWs were measured using the Automati c Routine for 
line E quivalent widths in stellar Spectra - ARES (ISousa et al.l 
l2007h . and they are also listed in Tabled 


The atomic lines of the elements Mg, V, Mn, and Cu have 
important HES that has to be taken i nto account. I n the case of 
those elements with lines not listed in lSteffenl(ll985h . we adopted 
the values of neighboring multiplets. An exception is the Mg line 
at 45785.29, for which no HSE information was available. Its gf 
value was obtained in the same way as those lines without HES 
listed in Table |2l This line is not strong (EW < 60 mA) and the 
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Fig. 3. Difference between the C abundances derived from the C 2 molecular bands (/15128 and /15165) and the C atomic line 
(/15380.34) as a function of the effective temperature (top panels) and of the S/N (bottom panels) for our sample of dwarfs (left 
panels), subgiants (middle panels), and giants (right panels). Stars with planets are represented by black symbols. 


error induced in the Mg abundance determination is of second 
order. The HFS of Ba, and also its isotopic splitting, are of some 
importance but only for the line at /16496 .91, being negligibl e 
for the lines at 45853.69 and 46141.73 (see lKorotin et alfcoi ih . 

Table |4] shows an excerpt of the derive photospheric param¬ 
eters for the 140 dwarf stars in our sample, together with their 
uncertainties. The complete tables, also for subdwarfs and gi¬ 
ants, are available in electronic form at the CDS. 


3.2. Abundance determination 

The abundances of the refractory elements (Na, Mg, Si, Ca, Ti, 
V, Mn, Ni, Cu, and Ba), with respect to Fe, were determined 
based on the equivalent widths of a selection of spectral lines 
chosen as being free from blends by inspection of t he Solar Flux 
Atlas and the solar line identifications catalogue of iMoore et al.l 

( 119 ^ . 

We used the abfind (in the case of elements without impor¬ 
tant HFS), and the blends (in the case of elements with HFS) 
drivers of the MOOG code to derive the abundances. These 
drivers require; (/) the model atmospheres, which were previ¬ 
ously computed using the photospheric parameters derived for 
each star; and (ii) a list of the selected lines, containing the 
respective wavelength (or the hyperfine structure), the atomic 
number, the lower-level excitation potential, the revised values 
of gf, and the measured equivalent widths. Again, using another 
routine for an automated calculation developed by us, we de¬ 
rived the abundances for each element in turn, excluding when 
needed (by means of a sigma-clipping procedure) spectral lines 
that provided a bad abundance determination. 

For the volatile elements (C, N, and O), and for Na as well, 
we used the spectral synthesis technique, also performed with 
the MOOG code, to reproduce the observed spectra in regions 
containing molecular lines of electronic-vibrational band heads 
of the C 2 Swan system (centered at 45128 and 45165), of the 
CN blue system (centered at 44215), and to atomic lines of C 
(45380.34), O (46300.30), and Na (46154.23 and 46160.75). 
Again, as in Paper 1, we preferred not to use the C atomic line at 


45052.17 because it is blended with a strong Fe line, which may 
affect the abundance determination. 

The synthetic spectra were computed in wavelength steps 
of 0.01 A, also considering the continuum opacity contribution 
in ranges of 0.5 A, line-broadening corrections (from velocity 
fields and instrumental broadening), and the limb darkening ef¬ 
fect. Some parameters of atomic and molecul ar spectral line s 
were taken from the VALD database and from iKuruczl (Il995h . 
respectively. Additional molecules that contribute to the spectral 
line formation in the studied wavelength regions are MgH (at 
45128, 45165, and for the C atomic line), CH (at 44215), and 
CN (for the Na and O atomic lines), though their contributions 
are relatively small. The gf values of some atomic and molecu¬ 
lar lines were revised to fit the solar spectrum, which was taken 
as a reference in our differential chemical analysis. The oxygen 
line at 46300.30 is blended with a Ni line at 46300.34 (logg/ = 
-2.3098). Therefore, we derived the O abundances by taking the 
Ni abundance of each star into account when fitting the observed 
profile. Moreover, to account for the presence of t elluric lines in 
this re gion, we adopted the telluric spectrum by [Wallace et al.l 
(I 2 OI ih . with the resolution degraded to match our spectra. 

To account for spectral line broadening from velocity fields, 
as in Paper 1, we adopted the same definition of Tbioad (also listed 
in Table |4]i, which is a composite of velocity fields, such as ro¬ 
tation and macro-turbulent velocities. The instrumental broad¬ 
ening was estimated from the width of thorium lines, and the 
linear coefficient u of the stellar limb darkening was taken from 
iDfaz-CordovesI (1 1 9951) . 

Figure |2] shows the spectral synthesis applied to three of the 
regions used for abundance determination, the C 2 (45165) and 
CN (44215) molecular band heads and the O atomic line at 
46300.30, for two dwarf stars with different spectral S/N. The 
synthetic spectra were resampled in steps of 0.05 A to consis¬ 
tently match the wavelength scale of the observed spectra. For 
examples of spectral synthesis applied to other regions and to 
other stars, we refer the reader to Fig. 3 and 5 of Paper 1. 

For a comparison between the results provided by molecular- 
and atomic-line C abundances, Fig.|3]shows the differences plot¬ 
ted as a function of Teff and of S/N for our subsamples of dwarfs. 
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Fig. 4. Evolutionary tracks (magenta lines) and isochrones (blue lines) from Yonsei-Yale, calculated for different values of metallic- 
ity, showing how we derived the stellar masses and ages. The ranges of age (increasing from left to right, HD 94518 excepted) and 
mass (increasing from right to left, HD 94518 excepted) of the models plotted are indicated on the top of each panel. 


subgiants, and giants. Besides a slight underabundance (within 
errors, however) of molecular with respect to atomic C abun¬ 
dances for cool dwarfs, no systematic difference is observed in 
the whole range of or S/N. 

We estimated the uncertainties in the same way as in Paper 1, 
i.e., using our routine that includes the error propagation of the 
input parameters used by MOOG; the photospheric parameters 
Teff, logg, [Fe/H], and and the broadening velocity Vbroad- 
Each one in turn, the input parameters are iteratively changed by 
their errors, and new values of abundance are computed. The dif¬ 
ference between new and best determination gives the error due 
to each parameter uncertainty, and (r([X/Fe]) is a quadratic sum 
of individual contributions. The error in Ebroad was estimated to 
be of the order of 1 km s ' or smaller. 

Table|4]shows an excerpt with the derived [Ti/Fe] abundance 
ratios for the 140 dwarf stars. The complete tables with abun¬ 
dances of all the other elements, also for subdwarfs and giants, 
are available in electronic form at the CDS. The parameters of 
atomic and molecular lines are listed in Tables 5 and 6 , which 
are only available in electronic form at the CDS as well. Table 5 
contains the wavelength of the spectral feature, the atomic and 
molecular line identification, the lower-level excitation potential, 
the gf values, and the dissociation energy Do (for molecular fea¬ 
tures only). Table 6 (strong atomic lines) contains the same in¬ 
formation of Table 5, except the dissociation energy parameter. 

4. Evolutionary and kinematic parameters 

4.1. Stellar masses and ages 

We derived masses and ages through interpolations in t he grids 
of Yonsei-Yale (Y^) ev olutionary tracks and isochrones (lYi et alJ 
l200lt iKim et al]l2002h on the HR diagram. The determination 
of these parameters takes into account the metallicity of each 
star. The stars’ luminosities used in the diagrams were calculated 


with parallaxes taken fro m the new reduction of the Hipparcos 
data (Ivan Leeuwerj2007h . and bolometri c corrections ( EC) from 
iFloweil (Il996h . which were revised by iTorresI (120101) . For the 
Sun, we adopt BCq - -0.07 and = 4.75. 

A few examples for different metallicities are illustrated 
in Fig. 0] To reproduce the Sun’s position in the HR 
diagrams, adopting - 5777 K and age = 4.53 Gy 
(iGuenther & Demarauel Il997h . the evolutionary tracks and 
isochrones were displaced in log(7’eff) by 0.001628 (~22 K in 
Teff) and in log(L/L 0 ) by 0.011. By doing so, we derived for 
the Sun a mass M - 1.003 + 0.004 Mq, and an age -4.31 + 
0.17 Gyr (see Fig.|4^). The figure also shows; a metal-poor and 
young star (HD 106516, Fig. HJr); a metal-rich one (HD 12661, 
Fig Et); a young star that is outside the lower border of the mod¬ 
els (HD 26767, Fig.|4}l); a low-mass and old star (HD 208313, 
Fig. 0^); and another star outside the limits of the models, this 
time in the region of low-mass and old stars (HD 94518, Fig.0f). 

The star HD 26767 is known to be a member of the Hyades 
cluster, of 625 Myr, and has been studied in several works in 
the past decades. Other stars either very young or outside the 
lower border of the models have been identified as members of 
the Hyades, the Ursa Major (~500 Myr), or the Pleiades (20- 
150 Myr) moving groups, or even classified as possible sub¬ 
dwarfs. In the case of old stars outside the model limits, they 
are likely low-mass stars with age no ol der than the age of the 
Universe (13.8 Gvr. lKomatsu et aril201 ih . 

Table Q lists the derived masses and ages, together with the 
parameters used to compute the luminosities. The uncertainties 
were estimated taking into account the error propagation of the 
photospheric parameters and of the luminosity, using cr(7r) from 
Hipparcos, and adopting cr{V) - 0.01 and cr{BC) - 0.005. In 
view of the asymmetric distribution of the isochrones on the HR 
diagram, the table shows the lower and the upper uncertainties 
on age. 
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Fig. 5. Abundance ratios as a function of metallicity for our sample of 120 dwarf stars of the thin disc population. Stars with 
planets are represented by black symbols. For sodium, the abundance ratios from both spectral synthesis and equivalent widths 
determination are plotted. The 95% confidence intervals (solid lines) around linear regressions, and their respective correlation 
coefficients (r) are also shown. The linear regressions were derived either for [Fe/H] < 0, [Fe/H] > 0, or for the whole metallicity 
range, depending on which one has the more significant slope. 


4.2. Kinematic properties 


In the same way as in Paper 1, we grouped our stars accord¬ 
ing to their main population membership by computing the 
probability that a given star belongs to the Galactic thin disc, 
thick di sc, halo, or to transition a pop ulation. First, using equa¬ 
tions of iJohnson & Soderblotij (Il987h . we computed the space- 
velocity components I/lsr. Flsr. and VkLSRwith respect to the 
Local Standard of Rest (LSR) using parallaxes and proper mo¬ 
tions from the new reduction of Hipparcos, and radial velocities 
from the ELOD I E arch ive, from the Simbad database, and from 
iHolmberg et al.l (l2007h . For the Sun, we adopted I/t .sr = 10.0 , 
Flsr = 5.3, and Wlsr = 7.2 km s~* (Dehnen & Binnevlll998h . 
Then, employing the equations of iReddv et alJ ( 20061) . we com¬ 
puted the probability that a star belongs to one of the three pop¬ 
ulations. A probability Pthin. Pthick, or Phaio greater than or equal 
to 70% classifies the star as a thin, thick, or an halo member, re¬ 
spectively. If the probability is in-between, the star is classified 
as belonging either to the thin/thick or to the thick/halo transition 
group. 


Our sample of 309 stars contains 281 thin-disc (out of which 
120 dwarfs, 25 subgiants, and 136 giants), and 16 thick-disc (10 
dwarfs, 2 subgiants, and 4 giants) population members. None of 
them belongs to the halo, and 13(10 dwarfs, 2 subgiants, and 1 
giant) belong to a thin/thick transition population. The member¬ 
ship of each star is indicated in Table |7] 


5. Results and discussion 

5.1. Abundance trends with metaiiicity 

The [X/Fe] abundance ratios as a function of metallicity are plot¬ 
ted in Figs.|5]|6l and|7]for dwarfs, subgiants, and giants, respec¬ 
tively. Since most of our sample stars belong to the thin disc, for 
homogeneity, only these stars were used in the abundance anal¬ 
ysis. They are all listed in the tables, but only thin-disc members 
are displayed in the figures. Besides iron, titanium is the only 
element in our study having spectral lines of both neutral and 
singly ionized atoms. We have found no significant difference 
between Ti i and Ti ii, hence we adopted an average value of all 
lines together. 
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Fig. 6. The same as Fig. |5]but showing the abundance ratios of 25 subgiant stars of the thin disc population. 


Despite the small number of metal-poor ([Fe/H] < -0.3 dex) 
or metal-rich ([Fe/H] > 0.3 dex) stars in this sample, the trends 
of the abundance ratios with metallicity are in good agreement 
with those published in the literature concerning the nucleosyn- 
thetic origin of the elements and the evolution of their abun¬ 
dances w ith time. Regard i ng the abun dance trends for d warfs 
see, e.g.. lAdibekvan et ^ (l201 2b) andikanget al.l (1201 ih . For 
giant stars see. e.g.. lTakeda et al. ( 20081) . The abundance trends 
for subgiants are very similar to those for dwarfs for all the in¬ 
vestigated elements. A few cases concerning the abundances of 
the three subsamples deserve particular attention, and we discuss 
them in the following paragraphs. 


There are systematic differences in the abundance ratios of 
carbon, nitrogen, and sodium when the subsamples of dwarfs, 
subgiants, and giants are compared with each other; in the range 
of solar metallicity (from -0.1 to 0.1 dex) giant stars are under- 
abundant in [C/Fe] by -0.17 + 0.07 dex, and overabundant in 
[N/Fe] by 0.31 + 0.10 dex and in [Na/Fe] by 0.16 + 0.08 dex. 
The abundance trends for these elements confirm our previous 
result f or carbon discussed in Paper 1, and agree with the re¬ 
sults of iTakeda et al.l (l2008h for sodium and carbon. Concerning 
oxygen, we found [O/Fe] = 0.07 + 0.07 dex in the solar metal¬ 
licity ran ge. Therefore , we do not confirm the oxygen deficiency 
found bv lTakeda et al.l (l2008h . who used the forbidden [01] line 
at 5577 A. Their oxygen abundances were in any case proba¬ 


bly underestimated as affirmed by th e authors. Our results also 
agree with those of iLiu et^ (1201 Ol) . who found that [C/Fe] is 
depleted by 0.13 dex, [Na/Fe] is overabundant by 0.1 dex, and 
[O/Fe] was not altered after the first dredge-up. From a theoreti¬ 
cal point of view, changes in the surface abundances of evolved 
stars is expected for C and N, an d unexpected for O according to 
quantitative predictions made bv llbenl([l99lh . 

A commonly used explanation for such systematic differ¬ 
ences is based on the effect of mixing processes that dredge up 
C-poor and N- and Na-rich material (produced by CN and NeNa 
cycles) to the surface of evolved stars. On the other hand, the 
use of spectral lines contaminated by blends is also a possibility. 
In order to test this, we derived the sodium abundances using 
the two techniques described in Sect. 13.21 from the equivalent 
widths and from spectral synthesis. The available Na lines in our 
spectral coverage are affected by C 2 and CN molecular bands. 
Performing spectral synthesis allows the inclusion of the already 
derived C and N abundances and, therefore, to account for the 
presence of such bands. However, doing so, the [Na/Fe] over¬ 
abundance in giants is only slightly reduced, not enough to can¬ 
cel out the larger systematic difference with respect to dwarfs, 
favoring the mixing-process hypothesis. Another result that fa¬ 
vors this hypothesis is that [C/Fe], [N/Fe], and [Na/Fe] have a 
trend with mass: a negative trend with increasing mass in the 
case of carbon, and positive ones for nitrogen and sodium (see 
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Fig. 7. The same as Fig. |5]but showing the abundance ratios of 136 giant stars of the thin disc population. The stars indicated are 
discussed in Sect lS.l.fl 


Fig. iJ. More massive stars have evolved faster and, therefore, 
have become C-poorer and N- and Na-richer than stars with 
masses close to the solar value. The oxygen abundances remain 
constant and independent of mass, supporting the above discus¬ 
sion. 

The nitrogen abundances of our giant stars show a clear 
dependence on the effective temperature (see Fig. |9] top right 
panel). One possibility is that this dependence is related to the 
continuum definition when using the CN molecular band. In this 
case, the reader would ask if the trend in [N/Fe] vs. [Fe/H] is 
also related to this dependence on Teff, due to a bad contin¬ 
uum definition. However, if we limit the effective temperature 
within the almost flat region seen in the [N/Fe] vs. Teff plane, 
from about 4600 to 5000 K, the trend of [N/Fe] with metallicity 
still remains. Moreover, the effective temperature of our sample 
of giants increases with increasing mass, i.e., it seems that more 
evolved stars have higher temperatures. Therefore, we believe 
that the [N/Fe] dependence on is probably related to evolu¬ 
tion effects instead of to bad continuum definition. 

Concerning the car bon abundances in dwarf sta rs, our result 
is in line with those of iDelgado Mena et al.l (1201 Oh in the range 
of solar metallicities. They found, however, a steeper trend in 
[C/Fe] vs. [Fe/H] in the range of metal-poor stars, and prac¬ 
tically a flat distribution of metal-rich dwarfs. On the other 


hand, [Gonzalez Hernandez et al.l (1201 Ol 1201 3bh derived [C/Fe] 
ratios continuously decreasing with increasing [Fe/H] in the 
whole range of metallicities for their sample of dwarfs and solar 
analogs. 

For the nitrogen abundances of dwarf stars, we derived a 
slope of 0.31 + 0.06 in the [N/Fe] vs. [Fe/H] plane for dwarfs 
with and without p lanets fitted together. T his is quite steeper than 
the one derived bv lEcuvillon et al.l (l2004h for their sample of so¬ 
lar type stars with and without planets (0.10 + 0.05) based on 
the NH band at 3360 A. We checked whether our [N/Fe] cor¬ 
relates with the effective temperature or the surface gravity, and 
no significant trend is observed in both cases. However, it seems 
that our dwarf stars with temperature higher than the solar value 
either have [N/Fe] close to zero or are slightly overabundant, 
whereas those cooler than the Sun are slightly underabundant 
(see Fig. |9] top left panel). By performing the same fit in the 
[N/Fe] vs. [Fe/H] plane but only to stars with Teff < 5800 K we 
derive a slope of 0.28 + 0.06, a value that is less steeper but still 
significantly high. 


The decreasing [O/Fe] with increasing metallicity observed 
for our sample of dwarfs agrees with the trends derived by 
[Delgado Mena et ^ ([201 Oh and by [Gonzalez Hernandez et al.l 
( 2010 [ 2013bh . and also with the recent results by 
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Fig. 8. Abundance ratios of C, N, O, and Na as a function of mass for our sample of dwarfs (left panels), subgiants (middle panels), 
and giants (right panels). Stars with planets are represented by black symbols. 


iBertran de Lis et ^ (1201 5h for the O abundances based on 
the Oi 26158.17 and [Oi] 26300.30 atomic lines. 

The distributions of [X/H] are shifted towards higher abun¬ 
dances (see Fig. [O) due to the metal-rich nature of stars host¬ 
ing giant plan ets as already known and discussed in the litera- 
ture (see, e.g.. iNeves et al.fl2009l: lAdibekvan et al.l 12012^ . The 
exceptions in the current study are oxygen and barium, which 
[X/H] abundances show no clear enrichment in giant planet hosts 
(see discussion in Sect. |6]l. In other words, most of these el¬ 
ements follows the same behavior of iron and are, for some 
reason, also linked to the planetary formation process. It is 
worth noting that only stars with giant planets (mass > 50 Earth 
masses) were considered and compared with only stars for which 
no giant planet has been detected. We conducted an exhaustive 
search in the literature in order to classify each star and create 
these two subsamples. They are identified in the remarks column 
of Table|7] together with the respective reference. In addition, we 
note that we performed this comparison only for dwarf stars of 
the thin disc population. Evolved stars are normally not included 
in the surveys of planet search, because they are well known for 
their intrinsic radial-velocity variation due to pulsation and/or 
jitter (distortion in the observed signal caused by stellar mag¬ 
netic activity). Subgiant stars could be added to the subsample 
of dwarfs, but we preferred to keep the homogeneity. 

In order to quantify the observed trends with metallicity, 
we fitted linear regressions in the range of stars metal-poorer 
or metal-richer than the Sun, or in the whole range of metallic¬ 
ity, depending on which one provides a higher correlation. The 
abundance ratios were then subtracted by the predictions of these 
linear fits trying to remove (or at least minimize) the effects of 
the chemical evolution of the Galaxy’s thin disc. The abundances 


corrected from such effects are used and discussed in Sect. 15.31 
where we derive the trends with condensation temperature of the 
elements. 


5.1.1. Chemically peculiar stars 

Most of the dwarf and subgiant stars plotted in Eigs. |5] and |6] 
follow the abundance trends with metallicity. A particular ex¬ 
ception seems to be the vanadium abundance ratios for some 
of them. Such stars have [V/Ee] higher than what is expected 
for their metallicity. By plotting the abundance ratios of the el¬ 
ements as a function of the other stellar parameters, we noticed 
that there is a systematic positive trend of [V/Ee] as a function 
of the effective temperature for Teff > 6000 K for our sample 
of dwarfs (see Eig. |9] bottom left panel). Eor this reason, in the 
following sections we discuss the results obtained by including 
or excluding the vanadium abundances of stars with tempera¬ 
ture higher than this limit. No clear trend with temperature for 
the other elements is observed, and of the abundance ratios with 
the other stellar parameters (except age) neither. This is the case 
at least for dwarfs and subgiants, because the elemental abun¬ 
dances in giants may have been affected by stellar evolution (see 
discussion in Sect. lO . 

Concerning the chemical peculiarity of giant stars, two of 
our sample stars deserve special given their position in the 
[Ba/Ee] vs. [Ee/H] diagram of Eig. |7] Both HD 104979 and 
HD 202109 appear enriche d in barium comp ared to other giants 
with the same metallicity. IWilliamsI (1197 Ih derived a [Ba/Ee] 
= 1.1 dex for HD 104979 and first suggested it to be a Ban 
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Fig. 9. [N/Fe] and [V/Fe] as a function of the effective temperature for our sample of dwarfs (left panels), subgiants (middle panels), 
and giants (right panels). Stars with planets are represented by black symbols. 
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Fig. 10. Abundance ratios as a function of age for our subsample of 120 dwarf stars of the thin disc population. Stars with planets 
are indicated by black symbols. For sodium, the abundance ratios from both equivalent widths and spectral synthesis determination 
are plotted. The vertical dashed line indicates the adopted solar age (4.53 Gyr). The 95% confidence intervals (solid lines) around 
linear regressions, and their respective correlation coefficients (r) are also shown. The linear regressions were derived either for age 
< 4.53 Gyr, age > 4.53 Gyr, or for the whole range in age, depending on which one has the more significant slope. 
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Fig. 11. Abundance differences [X/Fe]sun - [X/Fe]stars as a function of the condensation temperature of the elements for some 
of our dwarf stars of the thin disc population. We show the trends before (left panels) and after applying the corrections for the 
Galactic chemical evolution effects based on age (middle panels) and on metallicity (right panels). The standard deviations of the 
linear regressions (solid lines) for Tc > 900 K and their respective slopes are also shown. The planet host star is indicated by black 
symbols. 


sta]0. More recent works have confirmed the barium over abun¬ 
dance for this star, namely [Ba/Fe] = 0.93 + 0.10 dex (Zac^ 
11994 and [Ba/Fe] = 0.54 + 0.31 dex (ISmilianic et al.ll2007h . 
both in agreement with our determination of [Ba/Fe] = 0.77 + 
0.06 dex. The barium sta r natur e of HD 202109 was first sug¬ 
gested by IChromev et al.l (1 19691) . who measured [Ba/Fe] = 0.8 
+ 0.3 dex, and afterward classified as mild ( marginal enri ched) 
barium stars: [Ba/Fe] = 0.41 ± 0.13 dex by Zaci (1^994 . and 
[Ba/Fe] = 0.31 ± 0.26 dex by ISmilianic et alJ (l2007l) .We de¬ 
rived [Ba/Fe] = 0.40 ± 0.09 dex, a result that agree with the mild 
barium star classification. 


5.2. Abundance trends with age 

Figure [10] shows the [X/Fe] abundance ratios as a function of 
age for the subsample of dwarf stars of the thin disc popula¬ 
tion. Similarly to what is described in Sect. IS.ll in the case of the 
metallicity, we also quantified the observed trends with age by 
fitting linear regressions in the range of stars younger or older 
than the Sun, or in the whole range of age, depending on which 
one provides a higher correlation. Again, the abundance ratios 
were subtracted by the predictions of these linear fits in order to 
remove (or at least minimize), using another approach, the ef¬ 
fects of the chemical evolution of the Galaxy. The abundances 
corrected from such effects are used and discussed in Sect. 15.31 

^ A Ban star originates in the following scenario: when the more 
massive component of a binary system evolves as a thermally pulsing 
asymptotic giant branch star, the material produced in the He-buming 
envelope, enriched in s-process elements, is dredged up to the surface 
and then accreted by the companion by wind mass transfer. The initially 
more massive star finishes as a white dwarf whereas the companion 
becomes a primary barium star. 


where we derive the trends with condensation temperature of the 
elements. 

Depending on the position of the star in the HR diagram, the 
errors in age may be quite high, and the age determination is 
thus less reliable. By means of the distribution of the uncertain¬ 
ties on this parameter, we set a 2cr confidence limit up to which 
about 95% of the stars are found to have errors smaller than that 
limit, namely ~4 Gyr. Hence, only stars with (T(age) < 4 Gyr 
was considered in the analysis. 

5.3. Abundance trends with condensation temperature 

During the last few years, several works have investigated the 
relation between the stellar chemical abundances and the con¬ 
densation temperature (Tc) of the elements in stars with and 
without planets. Aiming to minimize systematic errors, the abun¬ 
dance difference between the Sun and other stars, A[X/Fe], has 
been used, where the solar abundances are derived from, e.g., a 
Moon’s spectrum. 

On the one hand, it has been suggested that the observed 
trends of A[X/Fe] as a function of Tc may have some relation 
with the planetary formation, in particular, with the formation 
of ter restrial planets (see the review performed bv iRamirez et al.l 
l201(]h . On the other hand, other authors argued that this kind of 
relation practically vanis hes after removing the Galactic ch emi- 
cal evolution effects (see Gonzalez Hernandez et al.ll2013b[ and 
references therein). Indeed, positive slopes in these diagrams 
may be the result of having higher abundances of less refractory 
elements, such as Na and Cu, with respect to more refractory 
ones, such as Ti and Ca. In Fig.|5] both [Ti/Fe] and [Ca/Fe] have 
a slight decrease with increasing metallicity whereas [Na/Fe] 
and [Cu/Fe] clearly increase for higher metallicities. 
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Fig. 12. Slope of the regressions computed in the diagrams A[X/Fe] vs. Tc (Fig. [TTI) plotted as a function of the derived stellar 
parameters. We show possible trends before (left panels) and after applying the corrections for the Galactic chemical evolution 
effects based on age (middle panels) and on metallicity (right panels). Dwarf stars of the thin disc population are separated in solar 
analogs (red symbols) and non-solar analogs (black symbols). Stars with planets are indicated by crosses. The slopes of linear 
regressions (solid lines) applied to the solar analogs and their respective correlation coefficients (r) are also shown. 


In particular, some of the studies discussed in the papers 
above consider only solar twins (T^ff = 5777 + 100 K, [Fe/H] 
= 0.0 + 0.1 dex, log g - 4.44 ± 0.10 dex) or solar analogs (Teff = 
5777 + 200 K, [Fe/H] = 0.0 ± 0.2 dex, logjg ^ 4.44 ± 0.20 dex) 
in their analysis. Recently, lAdibekvan et al.l (l2014li have pub¬ 
lished the results of their study in which they searched for corre¬ 
lations of the A[X/Fe] vs. Tq trends with some stellar and orbital 
parameters. After correcting the abundances from the Galactic 
chemical evolution effects, they found that the observed trends 
correlate with age and anticorrelate with the surface gravity, and 
possibly with the mean Galactocentric distance, and that the cor¬ 
relations (or anticorrelations) are even steeper and more evident 
if only solar analogs are plotted. 

Here we also investigate the relations involving A[X/Fe] as 
a function of the condensation temperature of the elements (to 
compute [X/Fe]sun. we used the solar abundances derived from 
the spectrum of the sunlight reflected by the Moon). First, we 
fitted linear regressions for elements with Tc > 900 K, a limit 
normally used to dist inguish refractory fro m volatiles. Values of 
Tc were taken from iLodders et aP (l2009l) . In order to account 
for any effects that the Galactic chemical evolution may have on 


the observed trends, we applied the corrections based on both 
metallicity and age, as described previously. Figure [TT] shows a 
few examples. Table 8, available in electronic form at the CDS, 
lists the slopes of the linear regressions and the respective stan¬ 
dard deviations for all dwarfs. Only dwarf stars were used in 
this analysis because their abundances are normally not affected 
by stellar evolution. The table contains the slopes and their un¬ 
certainties, the correlation coefficients, and the standard devia¬ 
tions of the regressions before and after applying the corrections 
for the Galactic chemical evolution effects based on age and on 
metallicity. 

We then plot (Fig. fTSI) the respective slopes as a function of 
the stellar parameters T^s, [Fe/H], log g, mass, and age. In this 
figure, solar analogs are distinguished from non-solar analogs, 
and stars with planets from those without planets. Regarding the 
sample of solar analogs, the figure shows: i) a marginal corre¬ 
lation with metallicity (in the sense that metal-richer stars have 
more positive slopes) and with mass after removing the effects 
of chemical evolution based on age; however, these correlations 
completely disappear if the correction is based on metallicity; 
ii) a clear anticorrelation with the surface gravity in both sce- 
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narios of Galactic chemical evolution correction; Hi) a marginal 
correlation with age, slightly steeper after applying the correc¬ 
tion based on metallicity. T he directions of these co rrelations 
agree with those reported bv lAdibekvan et all (l2014ll . We can¬ 
not check the Tq trends with Galactocentric distance because 
we only have solar neighborhood stars (distance to the Sun < 
100 pc). Concerning the relation with the presence of planets, 
the small number of planet hosts in our sample of solar analogs 
prevents us from any conclusions. 

By excluding from the analysis the vanadium abundances of 
stars with > 6000 K, the correlations and anticorrelations 
obtained above are still the same, specially because this limit is 
already imposed when only solar analogs are used. 

6. Statistical analysis 

We ha ve used a number of statistical methods (see lBreiman et alJ 
1 19841) to explore the multivariate properties of our sample. We 
are interested in uncovering significant differences in the stel¬ 
lar properties for stars with and without a known giant planet. 
Hereafter we will use the indices 1 and 0 to distinguish these two 
samples, respectively. Moreover, the sample of dwarf stars will 
be prefixed by the letter “d”. Thus, samples dl and dO stand for 
dwarfs with and without a known giant planet, respectively. The 
analysis and the following discussion is restricted to the dwarfs 
because their sample is thought to be less affected by an obser¬ 
vation bias. 

6.1. Abundance distributions 

We compare the abundance distributions between the samples 
dO (red) and dl (blue) in Fig. [13] The boxes and whiskers sum¬ 
marize graphically percentiles of these abundance distributions. 
Particularly, the middle line inside each box marks the median 
for that distribution, whereas the upper and lower edge of the 
box indicate the interquartile range. We applied a 2cr clipping in 
the abundance data of each single element in order to get rid of 
outliers in the abundance space. Panel (a) shows the comparison 
for all dwarf stars after applying the sigma clipping, while panel 
(b) restricts the comparison to those stars having -0.25 < [Fe/H] 
< -1-0.35 (a range in common between the subsamples of stars 
with and without planets). Both comparisons depict that dwarf 
stars with giants planets (dl) have generally larger abundances 
in nearly all elements compared to dwarfs without giant planets 
(dO). We note that vanadium seems to be exceptionally abundant 
in dl stars, whereas no significant difference is seen between 
dl and dO for the oxygen and barium abundances. Concerning 
vanadium, the overabundance of dl with respect to dO does not 
change if only stars with Teff < 6000 K are used. 

We used the Welsch’s t-test to quantitatively assess whether 
the means of these distributions are significantly different from 
one another. Table [1] gives the test statistics for stars having 
-0.25 < [Fe/H] < -1-0.35, after applying 2cr clipping. Some p- 
values are very low (< 0.001), reinforcing the visual comparison 
provided by Fig.[T3 The mean of the abundance distribution of 
all elements, but O and Ba, are higher in dl stars, compared to 
dO stars. Oxygen abundances are not very well represented in 
our sample, and only 32 dO and 4 dl stars have measured O 
abundance; this could be one of the reasons for why the t-test 
gave a large p-value for the O abundance distribution compari¬ 
son. Notwithstanding, it is highly interesting, although not clear, 
why Ba is not particularly more abundant in dl stars. 

The t-tests described above can only be applied to the abun¬ 
dance distribution of a single element at a time. They can pin¬ 


point elements that are particularly more or less abundant in one 
of the tested distributions. To use the information encoded in the 
whole abundance space, we also applied a multivariate analysis 
of variance (MANOVA) to assess for the significance of differ¬ 
ences in the multivariate means (that is, the location of samples 
0 and 1 in the abundance space). According to this MANOVA, 
we can reject the null hypothesis that the dO and dl samples have 
similar locations in the abundance space {p - 0.006). 

6.2. Abundance difference per metaiiicity bin 

We have also considered a different approach to disentangle 
the effects of the chemical evolution of the Galaxy from the 
metallicity-giant planet connection. We have defined a crude 
abundance quantity iZ = [Mg/H] + [Si/H] + [Ca/H] + [Ti/H] H- 
[V/H] -H [Mn/H] -I- [Ni/H] -H [Cu/H] -i- [Ba/H]. Z does not have 
an exact physical meaning, on account of the sum of unsealed 
log abundances, but it can serve as enrichment index for the pur¬ 
pose of this analysis. For instance, it is clear that Z and [Fe/H] 
should be correlated as a result of the Galactic chemical evolu¬ 
tion. This index could only be defined for those stars having a 
complete set of bundances of the elements Mg to Ba, including 
Fe. CNO and Na abundances were not included in the Z defini¬ 
tion because this would lower substantially the final number of 
stars having a complete set of abundance measurements. 

The index was used in a linear model of the type 

Z=yeo+/8i[Fe/H] + e, (1) 

fitted by least squares to the dO stars. The residuals of the fit 
are given by e for each star. Following this, we calculated the 
residuals of the dl stars with respect to the same linear model 
fitted to the dO stars. A Welch’s t-test applied to these residual 
distributions indicate that the difference in their means are sig¬ 
nificant (ap-value p - 0.0023), with eji -fdo = 0.16 + 0.10. That 
is, for a given [Fe/H], dl stars have generally larger Z indices 
than dO stars. If we can consider the dO stars as a “null hypoth¬ 
esis” population, for which the elemental abundances only re¬ 
flect the Galactic chemical evolution, the residuals of Z for the 
dl stars using the linear model in Eq. [1] should have no chem¬ 
ical evolution bias. Again, if we calculate Z excluding [V/H], 
this result practically does not change (the new p-value would 
he p - 0.0026). On the other hand, if Z is calculated exclud¬ 
ing the [Ba/H] abundance ratios, we obtain p - 0.0005 and 
fdi - fdo = 0.22 + 0.12, i.e., the difference between the dO and dl 
populations is even larger. 

We defined Z by mixing elements with different nucleosyn- 
thetic origins, which might veil the real difference in abundance 
between stars with and without planets. However, the mean 
[X/H] distributions for a (Mg, Si, Ca, Ti) and for iron-peak (V, 
Cr, Mn, Co, Ni) elem ents in stars hosting giant p lanets are very 
similar, as reported bv iMata Sanchez et al.l (l2014t) . Indeed, if we 
calculate Z for a elements only we obtain eji -fdo = 0.08+0.06, 
whereas for our iron-peak elements (V, Mn, Ni) only we have 
fdi - fdo = 0.10 + 0.06. For stars without planets, the [X/H] dis¬ 
tributions are also very si milar, having a maximum a bundance at 
—0.1 for most species (iMata Sanchez et ^l2014l) . Therefore, 
we conclude that the mix of elements when calculating Z is not 
a bad approach, specially because the differences that we are 
comparing become more significant when the abundances of all 
these elements are added together. 


13 















R. da Silva, et al.: Abundances in dwarfs, subgiants and giants with and without planets 




Fig. 13. Boxplot for the comparison between the elemental abundances of dwarf stars with and without planets, show in blue and 
red, respectively. Two sigma clipping was applied to all abundance distributions before the comparison, to avoid abundance outliers. 
The top panel shows the distributions for dwarfs in the [Fe/H] abundance space, whereas the bottom panel restrict the comparison 
to those stars having -0.25 < [Fe/H] < +0.35. The original data points are shown beside each corresponding box to allow the 
visualization of the intrinsic scatter. 


full partition model without V 





Fig. 14. Classification trees for the planet-harboring status amongst dwarf stars in our sample. The nodes are annotated with a 
categorical binary variable that assumes values 1 or 0 for stars with and without a giant planet, respectively. The leaves are marked 
with statistics of the partitioned group: fraction of dl stars in the group (stars harboring a giant planets) and percentage of all dwarfs 
stars that were classified in that leaf. The top left tree is the full partition model for our sample. The other trees consider the case 
when two important elements are not known: V and Cu, one at a time or both. 
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6.3. Classification trees 

Assuming the general connection between chemical abundance 
and the presence of a giant planet, we can extend the reasoning 
in using the iron abundance as an indicator of the plausibility 
that the star harbors a planet. To do this, we applied a classifica¬ 
tion tree to a complete set of [X/H] abundances for the dwarfs 
in our sample. Classification trees give the optimal splitting of 
observations coming from different classes (or populations). In 
our problem, we consider two classes: stars with and without gi¬ 
ant planets. The tree is built through recursive partitioning of the 
abundance space into nested groups, aiming at terminal groups 
having highest purity (that is, having most, if not all, of its mem¬ 
bers coming from a single class). Each classification tree has 
nodes, branches and terminal leaves (or group). The nodes are 
annotated with some property (or variable) and a splitting rule 
which points to a branch or the other. The leaves are usually 
marked with statistics of the partitioned group. 

As previously explained, on account of the requirement to 
have a complete set of abundances, CNO and Na measurements 
were not used in this analysis. No sigma clipping to the abun¬ 
dances was used before applying the classification trees. We 
show in Fig. fTTlfour classification trees. The leaves are marked 
(and color-coded, for easy identification) with the most likely 
class of its members: 0 (red) if the star is not likely to harbor a 
giant planet, and 1 (blue) if the opposite holds. Two numbers are 
annotated bellow each terminal leaves, showing the fraction of 
dl stars in the group (stars harboring a giant planet) and the per¬ 
centage of all dwarfs stars that were classified in that leaf. Thus, 
the numbers 0.42 9% can be read as 9% of all dwarfs were clas¬ 
sified in that group and 42% of them harbors a giant planet. The 
upper left panel gives the classification tree when we consider 
all elements from Mg to Ba. Although this abundance space is 
formed by the measurements of 10 chemical species, only two 
of them are needed for an optimal partition: Cu and V. Copper 
abundances seems to be the best single indicator of the presence 
of a giant planet: only 8% of the stars having [Cu/H] < 0.15 dex 
harbors a giant planet. Vanadium comes second in refining this: 
78% of stars having [Cu/H] > 0.15 dex and [V/H] > 0.33 dex har¬ 
bors a giant planet. The three other classification trees in Fig. [14] 
shows the partition model when Cu, V or both variables are not 
used. If vanadium is taken out of the independent variable set, 
the whole 10-dimensional abundance space could still be parti¬ 
tioned using the Cu abundances alone. The purity is somewhat 
lower: 73% of the stars having 0.15 < [Cu/H] < 0.24 harbors 
a giant planet. If, on the other hand, we take Cu from the set 
of independent variables, V takes the lead as the most relevant 
element, followed by two more splittings using Ti and Ca. The 
two more meaningful (extreme) leaves give: i) 5% of stars hav¬ 
ing [V/H] < 0.19 dex harbors a giant planet; and ii) 71% of stars 
having simultaneously [V/H] > 0.19, [Ti/H] > 0.14 and [Ca/H] 
>0.15 dex harbors a giant planet. Taking both Cu and V out of 
the set of independent variables yields a classification tree based 
on Ti and Ca abundances that is an exact subtree from the third 
case (with no Cu abundances). 

6.4. Linear discriminant analysis 

A similar, but more sophisticated, approach can be done using a 
linear discriminant analysis (FDA), which finds the linear com¬ 
bination of observed properties (in our case, abundances) which 
best separates two or more classes of objects. The result is an 
index, or system of indices, that may act as a linear classifier. 
The difference from the classification tree is that this one uses 


Table 1. Probabilities of the t-test statistics and coefficients of 
the linear discriminant for the dO and dl classes. 


Abundance 

ratio 

A(d0) 

A(dl) 

p-value 

Linear 

discriminant 

[Fe/H] 

94 

23 

2.39 X 

-0.87 

[C/H] 

86 

21 

5.12 X 10-“ 


[N/H] 

51 

10 

6.82 X 10-3 


[0/H] 

36 

5 

9.21 X 10-‘ 


[Na/H] 

48 

10 

4.92 X 10-3 


[Mg/H] 

89 

20 

8.38 X 10-“ 

-5.43 

[Si/H] 

87 

23 

1.43 X 10-“ 

0.70 

[Ca/H] 

88 

20 

7.62 X IQ-'^ 

-0.14 

[Ti/V] 

79 

21 

2.76 X 10-“ 

8.23 

[V/H] 

85 

21 

2.03 X 10-’ 

2.27 

[Mn/H] 

91 

23 

3.56 X 10-3 

-3.50 

[Ni/H] 

88 

23 

1.64 X 10-“ 

5.62 

[Cu/H] 

81 

22 

6.76 X 10-“ 

0.99 

[Ba/H] 

82 

22 

1.91 X 10-‘ 

-2.01 


Notes. The second and third columns gives the number of stars in each 
group. The p-value gives the probability that the mean of both dO and 
dl distributions are statistically similar. 



-2-10123 
Linear discriminant index 


Fig. 15. Histograms of the linear discriminant indices defined in 
Tabled] The color code blue and red is used to indicate stars with 
(top panel) and without planets (bottom panel), respectively. A 
black line shows the kernel density function for these indices, 
and a white arrow marks the position of the mean in each distri¬ 
bution. 


a hierarchical recursive splitting of single properties, aiming at 
classes with highest purity possible, while the FDA searches for 
rotations of the original multi-parametric space that leads to the 
best discrimination of the classes in the entire space. 

We applied the FDA to the complete abundance set of Mg 
to Ba for dwarfs, looking for the linear combination of abun¬ 
dances that could better discriminate the dO from dl samples. 
The resulting linear discriminant is given in the last column 
of Table [1] and in Fig. [15] As explained, the FDA finds the 
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best linear combination of the original abundances that discrim¬ 
inate amongst the tested classes. The numbers in the table are 
the coefficients of this linear discriminant. An index formed by 
-0.87[Fe/H] -5.43[Mg/H] -t-...- h0.99[Cu/H] -2.01 [Ba/H] can 
thus be used to discriminate dO from dl stars. If our sample can 
be taken as unbiased and representative of the stellar abundances 
of stars with and without giant planets, we can use this discrim¬ 
inant index in a predictive way, when the giant planet-harboring 
status of the star is unknown a priori. As pointed before, the con¬ 
cept is similar to using a classihcation tree, although now we use 
all abundances to build a single optimizer classiher. The power 
of this classiher can be seen in Fig.[T5j where we show the his¬ 
togram of the discriminant index. 

We explored all possible combinations of abundance ratios 
involving the elements from Mg to Ba, using the Welch’s test to 
assess interesting pairs. After considering all possibilities, we 
have found that [V/Ca] seems to be a good abundance ratio 
for discriminating dwarfs with planets: the average difference 
in [V/Ca] between the dO and dl groups amounts to 0.07 dex, 
leading to the rejection of the null hypothesis with p - 0.005. 

Kernel probability density functions for the relevant abun¬ 
dance ratio that allows the discrimination of dwarfs with and 
without giant planets are shown in Fig. (Tb] The color code is 
similar to that used before: blue and red for stars with and with¬ 
out planets, respectively. The hgure shows the [V/Ca] amongst 
dwarf stars. Although there is considerable intersection between 
the red and blue curves, it is both remarkable and puzzling that 
this abundance ratio seems to discriminate stars according to 
their planet-harboring status. 


7. Conclusions 

The inclusion in the current study of more stars, with and with¬ 
out detected planets, and more elements aimed to improve the 
work presented in Paper 1. A summary of our new results and 
conclusions is as follows: 

i) the systematic differences in the abundance ratios of C, N, 
and Na among dwarfs, subgiants, and giants seem to con- 
hrm that mixing processes, together with cycles of nucle- 
osynthetic reactions, do modify the abundances in the sur¬ 
face of evolved stars; 

ii) such mixing processes seem, indeed, to alter the Na photo- 
spheric abundances in giant stars considering that, even tank¬ 
ing into account the influence of C 2 and CN molecular fea¬ 
tures, a Na overabundance is still observed; 

Hi) the slopes in the A[X/Fe] vs. Tc diagrams show a correlation 
with age and an anticorrelation with the surface gravity, even 
after accounting for the effects of the Galactic chemical evo¬ 
lution; these results appear only when the sample is restricted 
to solar-analog s t ars; in other words, as also pointed out by 
lAdibekvan et alJ (l2014ll . older and more evolved stars are, in 
same way, connected to a lower amount of more refractory 
elements available for the star formation; no signihcant cor¬ 
relation with other stellar parameters, and no relation with 
the presence of planets are observed; 
iv) our statistical analysis has provided the following results: 

a) the overabundance of other elements in dwarf stars host¬ 
ing giant planets conhrms the previous statement that not 
only iron is linked to the planetary formation process, but 
also C, N, Na, Mg, Si, Ca, Ti, V, Mn, Ni, and Cu; 

b) for O and Ba there are no signihcant differences between 
the samples of dwarfs with and without giant planets; oxygen 



Fig. 16. Kernel probability density functions for the relevant 
abundance ratio, [V/Ca], that allows the discrimination of dwarfs 
with and without giant planets. The color code blue and red 
is used to indicate stars with and without planets, respectively. 
These abundance ratio was found to discriminate dwarfs with 
respect to their planet-harboring status, in a similar role to that 
played by the linear discriminant indices. 

is not very well represented in our analysis, but the result for 
barium is intriguing and deserves further investigation; 

c) vanadium is the element having the most signihcant 
abundance difference between giant planet hosts and single 
dwarfs, and this difference is still signihcant if only stars with 
Teff < 6000 K are used; 

d) by combining together the elements from Mg to Cu ac¬ 
cording to Eq. [1] we found a signihcant overabundance, for 
the same metallicity, in giant planet hosts in comparison with 
our sample of single dwarfs; 

e) based on the classihcation tree analysis we have found 
that abundances of copper and calcium are special discrim¬ 
inators of dwarfs with and without giant planets (over the 
[X/H] scale); 

/) by means of a linear discriminant analysis, we have found 
that [V/Ca] seems to be a good abundance ratio for discrim¬ 
inating dwarfs with giant planets. 
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Table 2. Atomic line parameters of the elements used in the analysis and EW?, measured in the degraded spectrum of the Solar Flux 
Atlas. 


A [A] 

Id. 

A 

[eV] 

log^/ 

EW 

[mA] 

4 [A] 

Id. 

T 

[eV] 

logg/ 

EW 

[mA] 

4 [A] 

Id. 

T 

[eV] 

logg/ 

EW 

[mA] 

5247.06 

Fei 

0.09 

-5.022 

64.0 

6265.14 

Fei 

2.18 

-2.588 

87.4 

5426.24 

Til 

0.02 

-2.778 

11.3 

5322.05 

Fei 

2.28 

-2.878 

62.6 

6380.75 

Fei 

4.19 

-1.315 

53.6 

5471.20 

Til 

1.44 

-1.558 

7.7 

5501.48 

Fei 

0.96 

-3.105 

121.2 

6498.94 

Fei 

0.96 

-4.608 

47.4 

5490.15 

Til 

1.46 

-0.994 

22.1 

5522.45 

Fei 

4.21 

-1.397 

46.1 

6608.03 

Fei 

2.28 

-3.945 

18.1 

5648.57 

Til 

2.49 

-0.378 

10.8 

5543.94 

Fei 

4.22 

-1.055 

64.7 

6627.55 

Fei 

4.55 

-1.476 

28.5 

5679.94 

Til 

2.47 

-0.690 

5.9 

5546.51 

Fei 

4.37 

-1.098 

54.6 

6703.57 

Fei 

2.76 

-3.010 

37.8 

5739.46 

Til 

2.25 

-0.802 

7.3 

5560.22 

Fei 

4.43 

-1.080 

52.8 

6726.67 

Fei 

4.61 

-1.060 

47.4 

5866.45 

Til 

1.07 

-0.812 

50.2 

5587.58 

Fei 

4.14 

-1.738 

32.0 

6733.16 

Fei 

4.64 

-1.434 

26.9 

6064.63 

Til 

1.05 

-1.888 

9.3 

5618.64 

Fei 

4.21 

-1.310 

50.9 

6750.16 

Fei 

2.42 

-2.638 

74.7 

6091.18 

Til 

2.27 

-0.422 

15.8 

5619.60 

Fei 

4.39 

-1.422 

36.2 

6752.71 

Fei 

4.64 

-1.206 

37.8 

6126.22 

Til 

1.07 

-1.352 

24.7 

5633.95 

Fei 

4.99 

-0.379 

69.0 

5234.63 

Fell 

3.22 

-2.196 

88.1 

6258.10 

Til 

1.44 

-0.440 

52.0 

5635.83 

Fei 

4.26 

-1.598 

33.6 

5325.56 

Fell 

3.22 

-3.237 

40.3 

4568.35 

Tin 

1.22 

-2.796 

31.4 

5638.27 

Fei 

4.22 

-0.828 

79.2 

5414.07 

Fell 

3.22 

-3.574 

26.8 

4583.42 

Tin 

1.16 

-2.934 

28.1 

5641.44 

Fei 

4.26 

-0.997 

66.6 

5425.25 

Fell 

3.20 

-3.235 

41.5 

4657.21 

Till 

1.24 

-2.312 

52.9 

5649.99 

Fei 

5.10 

-0.818 

35.0 

5991.38 

Fell 

3.15 

-3.568 

30.2 

4798.54 

Till 

1.08 

-2.664 

44.4 

5651.47 

Fei 

4.47 

-1.733 

19.6 

6084.11 

Fell 

3.20 

-3.770 

21.1 

5211.54 

Tin 

2.59 

-1.484 

33.3 

5652.32 

Fei 

4.26 

-1.700 

28.8 

6149.25 

Fell 

3.89 

-2.719 

36.5 

5336.78 

Tin 

1.58 

-1.638 

72.0 

5653.87 

Fei 

4.39 

-1.388 

38.0 

6247.56 

Fell 

3.89 

-2.290 

55.4 

5381.02 

Till 

1.57 

-1.866 

61.8 

5661.35 

Fei 

4.28 

-1.813 

23.3 

6369.46 

Fell 

2.89 

-4.150 

18.6 

5418.76 

Till 

1.58 

-2.106 

50.0 

5662.52 

Fei 

4.18 

-0.628 

95.5 

6416.93 

Fell 

3.89 

-2.618 

41.1 

5657.44 

Vi 

1.06 

- 

5.8 

5667.52 

Fei 

4.18 

-1.299 

53.3 

6432.69 

Fell 

2.89 

-3.545 

42.6 

5668.36 

Vi 

1.08 

- 

5.3 

5679.03 

Fei 

4.65 

-0.774 

60.3 

6456.39 

Fell 

3.90 

-2.133 

62.9 

5670.85 

Vi 

1.08 

- 

19.0 

5701.55 

Fei 

2.56 

-2.218 

84.5 

6154.23 

Nai 

2.10 

-1.619 

38.0 

5727.66 

Vi 

1.05 

- 

8.0 

5731.77 

Fei 

4.26 

-1.150 

57.6 

6160.75 

Nai 

2.10 

-1.419 

52.5 

6090.22 

Vi 

1.08 

- 

35.2 

5741.85 

Fei 

4.26 

-1.624 

32.5 

4571.10 

Mg I 

0.00 

- 

111.6 

6135.37 

Vi 

1.05 

- 

11.5 

5752.04 

Fei 

4.55 

-0.890 

58.3 

4730.04 

Mg I 

4.34 

- 

67.0 

6150.15 

Vi 

0.30 

- 

11.5 

5775.08 

Fei 

4.22 

-1.139 

60.4 

5711.10 

Mg I 

4.34 

- 

109.0 

6199.19 

Vi 

0.29 

- 

13.6 

5793.92 

Fei 

4.22 

-1.620 

34.7 

5785.29 

Mg I 

5.11 

-1.890 

51.0 

6216.36 

Vi 

0.28 

- 

36.4 

5806.73 

Fei 

4.61 

-0.895 

54.8 

5684.48 

Si I 

4.95 

-1.638 

62.8 

6274.66 

Vi 

0.27 

- 

9.8 

5809.22 

Fei 

3.88 

-1.617 

51.2 

5690.43 

Si I 

4.93 

-1.800 

52.9 

6285.17 

Vi 

0.28 

- 

9.4 

5814.81 

Fei 

4.28 

-1.835 

22.6 

5701.11 

Si I 

4.93 

-1.994 

40.4 

4626.54 

Mni 

4.71 

- 

29.5 

5852.22 

Fei 

4.55 

-1.183 

41.5 

5793.08 

Si I 

4.93 

-1.918 

45.2 

4739.11 

Mni 

2.94 

- 

65.0 

5855.08 

Fei 

4.61 

-1.519 

23.0 

6145.02 

Si I 

5.61 

-1.420 

40.2 

5004.89 

Mni 

2.92 

- 

12.2 

5856.09 

Fei 

4.29 

-1.558 

34.5 

5581.98 

Cai 

2.52 

-0.780 

97.0 

5394.67 

Mni 

0.00 

- 

78.6 

5862.36 

Fei 

4.55 

-0.423 

89.6 

5590.13 

Cai 

2.52 

-0.829 

93.5 

5399.48 

Mni 

3.85 

- 

39.2 

5905.68 

Fei 

4.65 

-0.787 

59.4 

5867.57 

Cai 

2.93 

-1.566 

26.6 

5413.68 

Mni 

3.86 

- 

23.5 

5916.26 

Fei 

2.45 

-2.922 

55.0 

6161.29 

Cai 

2.52 

-1.362 

59.6 

5420.35 

Mni 

2.14 

- 

83.6 

5927.79 

Fei 

4.65 

-1.068 

43.2 

6163.75 

Cai 

2.52 

-1.295 

63.8 

5432.55 

Mni 

0.00 

- 

51.6 

5929.68 

Fei 

4.55 

-1.207 

40.5 

6166.44 

Cai 

2.52 

-1.188 

70.6 

5537.77 

Mni 

2.19 

- 

36.4 

5930.19 

Fei 

4.65 

-0.336 

90.8 

6169.04 

Cai 

2.52 

-0.842 

94.6 

6013.50 

Mni 

3.07 

- 

87.0 

5934.66 

Fei 

3.93 

-1.139 

76.3 

6169.56 

Cai 

2.52 

-0.621 

112.3 

6021.80 

Mni 

3.07 

- 

91.2 

5983.69 

Fei 

4.55 

-0.839 

61.6 

6455.60 

Cai 

2.52 

-1.415 

57.3 

5032.72 

Nil 

3.90 

-1.178 

24.9 

5984.82 

Fei 

4.73 

-0.346 

85.1 

6499.65 

Cai 

2.52 

-0.954 

87.5 

5094.41 

Nil 

3.83 

-1.052 

33.3 

6024.06 

Fei 

4.55 

-0.142 

113.8 

4518.02 

Til 

0.83 

-0.394 

74.2 

5220.30 

Nil 

3.74 

-1.246 

28.7 

6027.06 

Fei 

4.08 

-1.195 

65.1 

4548.77 

Til 

0.83 

-0.438 

72.5 

5392.33 

Nil 

4.15 

-1.328 

12.7 

6056.01 

Fei 

4.73 

-0.507 

73.5 

4617.25 

Til 

1.75 

0.288 

66.3 

5435.87 

Nil 

1.99 

-2.440 

52.3 

6065.49 

Fei 

2.61 

-1.687 

119.6 

4758.12 

Til 

2.25 

0.312 

44.5 

5452.86 

Nil 

3.84 

-1.482 

16.6 

6079.01 

Fei 

4.65 

-1.015 

46.4 

4759.27 

Til 

2.25 

0.340 

46.0 

5494.88 

Nil 

4.10 

-1.055 

23.0 

6082.72 

Fei 

2.22 

-3.572 

34.7 

4778.26 

Til 

2.24 

-0.380 

15.6 

5587.85 

Nil 

1.93 

-2.438 

55.8 

6094.38 

Fei 

4.65 

-1.553 

20.6 

4926.15 

Til 

0.82 

-2.240 

6.1 

5625.31 

Nil 

4.09 

-0.714 

39.3 

6096.67 

Fei 

3.98 

-1.788 

38.1 

5022.87 

Til 

0.83 

-0.452 

75.4 

5637.13 

Nil 

4.09 

-0.846 

32.8 

6151.62 

Fei 

2.18 

-3.302 

50.1 

5071.47 

Til 

1.46 

-0.738 

31.6 

6176.81 

Nil 

4.09 

-0.270 

65.3 

6157.73 

Fei 

4.08 

-1.244 

63.2 

5113.45 

Til 

1.44 

-0.846 

27.8 

6177.24 

Nil 

1.83 

-3.497 

15.8 

6165.36 

Fei 

4.14 

-1.488 

46.1 

5145.46 

Til 

1.46 

-0.626 

37.3 

6186.71 

Nil 

4.10 

-0.886 

31.5 

6180.21 

Fei 

2.73 

-2.572 

60.1 

5147.48 

Til 

0.00 

-2.154 

33.5 

6378.26 

Nil 

4.15 

-0.812 

33.2 

6188.00 

Fei 

3.94 

-1.629 

48.5 

5152.19 

Til 

0.02 

-2.088 

35.7 

5782.14 

Cui 

1.64 

- 

80.5 

6200.32 

Fei 

2.61 

-2.405 

74.7 

5192.97 

Til 

0.02 

-1.012 

87.9 

5853.69 

Ban 

0.60 

-0.888 

63.3 

6226.74 

Fei 

3.88 

-2.064 

29.7 

5211.21 

Til 

0.84 

-2.120 

7.9 

6141.73 

Ban 

0.70 

0.043 

114.3 

6229.24 

Fei 

2.84 

-2.860 

39.8 

5295.78 

Til 

1.07 

-1.670 

12.5 






6240.65 

Fei 

2.22 

-3.298 

48.7 

5219.70 

Til 

0.02 

-2.250 

28.6 







Notes. Lines with missing gf values represent the elements for which the hyperfine structure was taken into account (Mg, V, Mn, and Cu) and the 
detailed line splitting is shown in Table [3 
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Table 3. The gf values for lines with hyperfine structure computed based on the degraded spectrum of the Solar Flux Atlas. 


d[A] 

logg/ 

d[A] 

logg/ 

d[A] 

logg/ 

Mg 1:4571.10 

Vi : 6274.66 

Mni: 5420.35 

4571.078 

-6.554 

6274.640 

-2.150 

5420.277 

-2.311 

4571.087 

-6.594 

6274.658 

-2.150 

5420.301 

-2.233 

4571.096 

-5.694 

6274.676 

-2.150 

5420.334 

-3.093 

Mgi: 4730.04 

Vi : 6285.17 

5420.376 

-1.986 

4730.031 

-3.254 

6285.147 

-2.160 

5420.429 

-1.898 

4730.038 

-3.294 

6285.165 

-2.160 

Mni: 5432.55 

4730.046 

-2.394 

6285.183 

-2.160 

5432.512 

-4.348 

Mg 1: 5711.10 

Mni: 4626.54 

5432.540 

-4.434 

5711.074 

-2.716 

4626.464 

-0.778 

5432.565 

-4.544 

5711.083 

-2.756 

4626.504 

0.017 

5432.584 

-4.689 

5711.091 

-1.856 

4626.530 

-0.231 

5432.598 

-4.783 

Vi : 5657.44 

4626.565 

0.133 

Mni: 5537.77 

5657.418 

-1.569 

4626.573 

-0.327 

5537.691 

-2.802 

5657.436 

-1.569 

Mni: 4739.11 

5537.710 

-2.689 

5657.454 

-1.569 

4739.099 

-1.189 

5537.738 

-2.653 

Vi : 5668.36 

4739.113 

-1.330 

5537.764 

-2.689 

5668.344 

-1.590 

4739.126 

-1.485 

5537.802 

-2.325 

5668.362 

-1.590 

4739.145 

-1.042 

Mni: 6013.50 

5668.380 

-1.590 

4739.167 

-2.392 

6013.474 

-0.704 

Vi : 5670.85 

Mni: 5004.89 

6013.486 

-0.914 

5670.833 

-0.970 

5004.878 

-2.249 

6013.501 

-1.045 

5670.851 

-0.970 

5004.892 

-2.390 

6013.519 

-0.724 

5670.869 

-0.970 

5004.905 

-2.545 

6013.537 

-1.302 

Vi : 5727.66 

5004.924 

-2.102 

Mni: 6021.80 

5727.643 

-1.434 

5004.946 

-3.452 

6021.764 

-1.386 

5727.661 

-1.434 

Mni: 5394.67 

6021.780 

-1.237 

5727.679 

-1.434 

5394.617 

-4.028 

6021.797 

-0.406 

Vi: 6090.22 

5394.645 

-4.114 

6021.806 

-0.623 

6090.234 

-0.629 

5394.670 

-4.224 

6021.814 

-0.488 

6090.216 

-0.629 

5394.689 

-4.369 

Cm: 5782.14 

6090.198 

-0.629 

5394.703 

-4.463 

5782.032 

-3.554 

Vi: 6135.37 

Mni: 5399.48 

5782.042 

-3.857 

6135.352 

-1.280 

5399.435 

-0.904 

5782.054 

-3.156 

6135.370 

-1.280 

5399.446 

-1.094 

5782.064 

-3.207 

6135.388 

-1.280 

5399.479 

-0.995 

5782.073 

-3.511 

Vi : 6150.15 

5399.502 

-0.616 

5782.084 

-2.810 

6150.136 

-2.030 

5399.536 

-1.245 

5782.086 

-3.156 

6150.154 

-2.030 

Mni: 5413.68 

5782.098 

-3.156 

6150.172 

-2.030 

5413.613 

-1.826 

5782.113 

-2.810 

Vi: 6199.19 

5413.653 

-1.031 

5782.124 

-2.810 

6199.168 

-1.960 

5413.679 

-1.279 

5782.153 

-2.709 

6199.186 

-1.960 

5413.714 

-0.915 

5782.173 

-2.363 

6199.204 

-1.960 

5413.722 

-1.375 



Vi: 6216.36 





6216.340 

-1.416 





6216.358 

-1.416 





6216.376 

-1.416 
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Table 4. Excerpt from the list of 140 dwarf stars with the photospheric parameters and [Ti/Fe] abundance ratios. 


Star 

Spectral 

type 

Abroad 

[km s“^] 

Teff ± cr 
[K] 

logg ± a- 

^ ± tr 
[kms^*] 

[Fe/H] ± a- 

[Ti/Fe] ± a 

BD+290503 

G5 

0.00 

5180 ±50 

4.26 ± 0.20 

0.86 ±0.31 

0.11 ± 0.07 

-0.07 ±0.10 

HD 10145 

G5V 

0.00 

5593 ± 43 

4.25 ±0.19 

0.78 ± 0.14 

0.01 ± 0.07 

0.03 ± 0.07 

HD 10307 

G1.5V 

0.00 

5913 ±49 

4.29 ± 0.22 

1.06 ±0.09 

0.09 ± 0.06 

-0.02 ± 0.07 

HD 10476 

K1V 

0.00 

5141 ± 31 

4.29 ± 0.18 

0.40 ± 0.25 

-0.05 ± 0.05 

0.02 ± 0.08 

HD 106116 

G4V 

0.00 

5620 ± 36 

4.22 ± 0.17 

0.84 ±0.11 

0.16 ±0.05 

-0.01 ± 0.09 

HD 106516 

F9V 

0.00 

6401 ± 158 

4.85 ± 0.39 

1.59 ±0.66 

-0.47 ± 0.12 

0.29 ±0.17 

HD 10780 

KOV 

0.00 

5339 ± 26 

4.40 ±0.15 

0.71 ± 0.13 

0.00 ± 0.04 

0.02 ± 0.06 

HD 108954 

F9V 

0.00 

6194 ± 60 

4.60 ± 0.23 

1.15 ±0.13 

0.08 ± 0.08 

0.03 ± 0.09 

HD 109358 

GOV 

0.00 

5875 ± 30 

4.42 ± 0.14 

1.21 ± 0.08 

-0.22 ± 0.04 

0.03 ± 0.06 

HD 11007 

F8V 

0.00 

6032 ± 44 

4.16 ± 0.26 

1.39 ±0.09 

-0.14 ±0.06 

0.03 ± 0.09 


Notes. The broadening velocity Vbmad is also shown. The whole tables for this and other elements, also for subgiants and giants, are available in 
electronic form at the CDS. 


Table 7. The same as Table |4]but showing photometric and evolutionary parameters together with the population membership. 


Star 

B-V 

Mv 

dist. 

[pc] 

BC 

logL 

mass 

[Mq] 

age 

[Gyr] 

Population 

group 

Remarks 

Ref. 

HD 10145 

0.691 

4.84 

37.3 

-0.106 

0.005 ± 0.026 

0.94 ± 0.03 

10.1+2^18 

thin 



HD 10476 

0.836 

5.82 

7.5 

-0.239 

-0.331 ±0.006 

0.78 ± 0.04 

17.U'*-2_7.7 

thin 

a 


HD 106116 

0.701 

4.79 

34.9 

-0.100 

0.025 ± 0.015 

1.00 ±0.03 

7.U2.5 

thin/thick 



HD 106516 

0.470 

4.36 

22.4 

0.011 

0.150 ±0.016 

1.05 ±0.05 

0.2+'Vo,2 

thick 



HD 108954 

0.568 

4.53 

21.8 

-0.010 

0.092 ± 0.008 

- 

- 

thin 

b 


HD 128311 

0.973 

6.38 

16.6 

-0.126 

-0.603 ± 0.015 

- 

0.50 

thin 

c, d (U) 

1,2 

HD 124292 

0.733 

6.18 

22.1 

-0.143 

-0.516 ±0.011 

- 

- 

thin 

e 

3 

HD 206860 

0.587 

4.74 

17.9 

-0.020 

0.017 ± 0.008 

- 

0.10 

thin 

c, d (P) 

1,2 

HD 25825 

0.593 

4.50 

45.9 

-0.021 

0.108 ±0.052 

- 

0.65 

thin 

d(H) 

4 

HD 42807 

0.663 

5.17 

18.0 

-0.070 

-0.140 ±0.008 

— 

2.70 

thin 

d(W) 

5 


Notes. Not considered in the abundance analysis because the age uncertainties are larger than 4 Gyr (see Sect. 14.11 : ® Mass and age not derived 
because the star is located outside the evolutionary tracks and isochrones used; Star hosting at least one giant planet; Member of a stellar 
moving group (U: Ursa Major; P: Pleiades; H: Hyades; W: Wolf 630); Subdwarf candidate. 

References. http://exoplanet.eu; ‘^^ dMontes et all ( 1200 Ih ; ^^ dSandage & Fout^ 119861) ; <"^ dEggenl 1195011 ; ^^ dBubar & Kingl 120101) 
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